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ABSTRACT 

The problem of the threshold law for electron-atom impact 
ionization is reconsidered as an extrapolation of inelastic cross 
sections through the ionization threshold. The cross sections are 
evaluated from a distorted wave matrix element, the final state 
of which describes the scattering from the N excited state of 
the target atom. The actual calculation is carried for the e-H 
system, and a model is introduced in which the r^ 2 ^ repulsion is 
replaced by (r^+r 2 ) ^ • This model is shown to preserve the 
essential properties of the problem while at the same time reducing 
the dimensionability of the Schrodinger equation. Nevertheless, 
the scattering equation is still very complex. It is dominated 
by the optical potential which is expanded in terms of eigen- 
spectrum of QHQ. It is shown by actual calculation that the lower 
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eigenvalues of this spectrum descend below the relevant inelastic 
thresholds; it follows rigorously that the optical potential con- 
tains repulsive terms. Analytical solutions of the final state 
wave function are obtained with several approximations of the 
optical potential: (o) omission of the optical potential (1) 

inclusion of the lowest term and dominant pole term (2) a closure 
approximation which depends on an effective energy for each 
threshold energy E„. The threshold law in all these cases is 
obtained. In the closure approximation the law depends on the 
sign and N dependence of £ - p , however it cannot be 

excluded that the difference in an oscillating function of N. 

In that case the derivative of the yield curve is an oscillating 
(but non-negative) function of the available energy E. A form 
of such a threshold law is suggested. 
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I . INTRODUCTION 

1 2 

In previous papers ' we have begun to consider the impact 
ionization problem from a completely quantum mechanical point of 
view. The touchstone of our understanding of that problem is the 
threshold law, and it is to that specific problem that we return. 

The insight that we tried to gain was by a study of the 
doubly excited (i.e. auto-detaching) states of the electron-atom 
system associated with ever higher principal quantum numbers of 
the target atom. The actual extrapolation procedure that was 
used, however, was through a summation of inelastic cross sections 
to such higher states, in which the final state wave function was 
taken as being of the form as the doubly excited state which 
minimized the energy. 

As reasonable as this procedure would appear, it is at best 
speculative, because the doubly excited states actually enter 
the equation for the final state scattering functions as specific 
terms in the optical potential. For each scattering function there 
are an infinity of optical potential terms plus direct potentials, 
not to mention coupling terms between various excited states that 
must in principle be considered. In the light of this complexity 
it is naive to expect that the final state scattering function 
is simply of the form of the lowed lying doubly excited state. 

Thus we here consider the scattering problem itself. First 
we define a model which we believe contains all the essentials 
of the electron-hydrogen ionization problem and yet greatly 
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reduces the mathematical complexity: we replace the electron- 


electron repulsion 2_ (in rydberg units which we use throughout) 

r 12 

by 2/(r 1 +r 2 ): 

1 . 1 

r 12 r l + r 2 


( 1 . 1 ) 


and thereby reduce the S-wave Schrodinger equation to a two- 

3 

dimensional partial differential equation . As a result the ex- 
cited spectrum of target states contain only s-states and loses 
the J l degeneracy associated with the complete hydrogenic spectrum. 
Nevertheless the long range dipole potential which the scattering 
particle sees is retained- in the model. These and other character- 
istics will become clear as we go along. 

In Section II we consider the scattering problem starting 
from a general close coupling expansion. We show that because of 
the nature of the spectrum of Q N HQ N that for N large the optical 
potential starts to contain repulsive terms even when all the 
coupling is included. This is our most important rigorous obser- 
vation. We shall also argue (Section IJ) that for purposes of 
evaluating inelastic scattering matrix elements, we can neglect 
the coupling terms, i . e . in effect we are considering a distorted 
wave approximation and that is our most important approximation. 

The direct potential (Hpp) problem is considered in Section 
III. Here we can introduce some benign approximations which 
allow analytic solutions to be given, which are nevertheless 
essential for a cogent analysis of what happens in the limit 
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N+°°. Basically these are the zero energy solutions in a Coulomb 
and in a dipole potential . 

The optical potential is examined in Section IV. We consider 
three approximations: a lowest term approximation; an effective 
intermediate state (dominant pole) approximation; and an effective 
energy or closure approximation. In Section V the threshold law 
for these various approximations is worked out, and some dis- 
cussion of the results is given including comparison with other 

4 

recent approaches to the problem based on Wannier . 
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II . A MODEL OF THE ELECTRON-HYDROGEN INTERACTION AND THE 
SCATTERING PROBLEM 

We consider the Schrodinger equation (rydberg units throughout) 


W N 


ET 


N 


( 2 . 1 ) 


for model corresponding to (1.1). 

Xi~ 


The Hamiltonian is given by 


/f 


- j. 2~r, _i 

n ^ ' 


- <v. 


— 2L 

n 


r. 


j l ..... 

r f + r 2. 3 


and we expand the solution in two parts 


*N = p * + Q N *N 
N-N N ' 


(2.3) 


corresponding to open channels, 

A/ 


P Tp --"2_ a >( rj + 

-Hr, ^ ^ ' 


and closed channels 


for a total energy E where 


(2 S) 


(w T/' = (}_+[) (2. S') 

1 A' 'n j rj/^ 


E N iE<E N+1 


( 2 . 6 ) 


with E n being the energy of N*"* 1 excited state of hydrogen 


E N “n^ 


(2.7) 


[We consider for the present only singlet solutions giving rise to 
only the + sign in (2.4) and (2.5).] The functions U^r) are to 
be determined; the target states 4> n (r) are S eigenstates of the 


6 



of the hydrogen atom: 


<f> IT)= 7 R (r) 
» r n 


As is by now well known, an equation equivalent to the 
Schrodinger equation can be derived^ for the open channel wave 
function P '{'. T 

[P HP - E] P^=0, 

where the Q-part of the optical potential is given by 

" p HQ. 1 Q HP 

° Pt E^OTfO 

For use in Appendix C we define the Q-space Green's function in 
the above equation . 


G Q = 


E-Q HQ 


( 2 . 8 ) 


( 2 . 9 ) 


( 2 . 10 a) 


(2.10b) 


Explicit forms for P. and Q. can be given as simple generalizations 

L 

of the formulas for N=1 

Q “ Qg Q 2 / (2.11a) 

where 

N 

Q h = 1 - Z <f>_ (i) ><♦_ (i) r (2.11b) 

1 n=l n n ' 

and as usual 

P = 1 - Q . ( 2 . 12 ) 

However we shall not need them, as our functions will be constructed 
to be manifestly in P or Q space. 

In this and many other contexts it is most convenient to ex- 
pand the optical potential in terms of the eigenfunctions of Q .HQ 
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Q. PQ T $ Nv ^fov^Nv' ( ' 2 ' 

where the eigenfunctions $ Nv are understood to be in Q space: 

Q $ = $ (2 . 

w - Nv Nv 

Using (2.13) we obtain the spectral representation of the ( Q-part etf) 
th^.optical potential 

tCA = ( 2 . 

E - £ 

'uJV 

The expansion (2.15) is not only useful, but it manifests many 

/r<0 rr> 

features of interest. For example in scattering^low lying states 

the fact that the numerator of (2.15) is positive definite taken 

together with the fact that the lowest states of QHQ are just 

slightly below the next inelastic threshold (=> E- <^ v <o for E^E^) 

implies that the optical potential is negative definite (i.e. 

attractive) and this is the basis for lower bound principles for 

6 7 

the scattering phase shifts ' . 


However, this is a situation that only obtains for low- 
lying N as Table I shows. There we have computed 

$ > 


2 2 
— <$ H$ > = — + <(S 

N N N n 2 W 1+ r 2 Nv 


for 


$ N fn ^ r l ) ^N^ r 2^ 


(2 


(2 


and we see that 


4 


' 5,-1' E N-2' E N-u whenever 


In other words the lowest eigen value associated with higher 

J_ T_ 

N’ states can descend below the n threshold, and when this 
happens the contribution of those terms to the optical potential 


13) 

14) 

15) 


.16) 

.17) 
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is repulsive. This does not prove that the effect of the whole 
optical potential will be repulsive, but it does suggest that it 
may be repulsive, and that in any event its effect will have to 
be considered very carefully. 


On the right hand side of Table I we have given similar 
results for the full interaction, in which case 

refers to a configuration interaction wave function'*' 

^ C N£ ) (r l ) ^N£ (r 2 ) P £ ( cos8 i2> 
and the linear combination giving the lowest energy is given 
(j=l) . Details of this calculation are given in Ref. 1. The 
point of showing those results is to demonstrate that the ' lower- 
ing of QHQ eigenvalues below lower N states is a property of the 
full e-H problem and not simply of the model. Indeed the Table 
shows that the model is remarkably accurate. 


Finally it should be realized that in the model (and corres- 
pondingly in the complete interaction case) there are many other 
linear independent functions in Q space, for example 


2 


-'A- 


[^v tr,iR v/, rx) " J , 


fVj N f I 

which have similar type of spectral properties going over finally 
to the purely dipole type states (labelled ^ in Ref. 1) in 
which the outer electron sees the induced dipole moment from the 
inner electron and the nucleus. Here too, there are an infinity 
of states but that spectrum probably always remains between E„ 

and E. t 

N. 


J N-1 
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III. THE DIRECT POTENTIAL PROBLEM 


The Hpp problem, i.e., 

[P HP - E n ]PT n =o (3.1) 

is itself a complicated problem by virtue both of the coupling 
between different open channels as well as the exchange terms 
associated with The latter, however, involve the same type 

of integral terms as those coming from the optical potential 
without involving the small energy denominators . Thus they are 
negligible in this context (although it should be recalled that 
they are essential even for qualitative purposes in low energy 
elastic scattering from the ground state to give the right nodal 
structure to the scattered orbital) . 

The coupling terms in (3.1) involve terms of the form 

V nA r i> U m< r l> and assumi ng n and m are of the order of N then 

V n m ^ ^ for a11 values of r. In perturbation theory these 

N 

potentials are to be divided by the energy differences which 

are also of the order ~ (n-m) . But the energy differences are of 

N 3 

both signs, thus it is not unreasonable to assume that a kind of 
random phase phenomenon will ensue in which the various terms will 
have a cancelling effect on each other. Furthermore, it must be 
recalled that the physical distance between the various N shells, 
< N/r/N>-<N+l/r/N+l>«N, actually increases with N. Finally it 
should be realized that the wave function we are attempting to 
calculate is to be used in an integral expression for the in- 
elastic amplitude. This is consistent with the philosophy of the 
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distorted wave approximation that the integral expression corrects 
to some extent for the inadequacies of the approximations of the 
wave functions that one puts into it. None of these arguments, 
however, is intended to imply that the omission has been rigor- 
ously justified. 


The Hpp equation becomes in this approximation 

d^ - V + k ^kwr) = o 

— V N,N k N J N lr) ° 

dr 2 


[ 


where 


and 


k N = E ‘ E N 


V. 


N 




r l r 2 +r 


4v nn 

r l 


(r) 


Little v NN is then the diagonal element of the electron-electron 
repulsion and it alone survives in off-diagonal elements 


1/ (r} = or (r) s* <o? (rp) / -i- ) (p \ 

A )M T /V ' ' C+r ( 'm / I 

*• I 

Although the potentials in their entirety are complicated their 
effect in our application can be simply approximated by (N=M) 





r>r 0 


( 3 . 2 ) 

( 3 . 3 ) 

( 3 . 4 ) 

( 3 . 5 ) 

0 . 6 , 

(3-7) 
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(3.8) 


where r is the mean radius of the N state 
0 

r Q = <N|r|U> = | N 2 , 


and the dipole moment ; 


- 2 = + - 

K lr+f- 3 _ _ 2 

r 


~> 


b N + 0 (r -3 ) 


whe r-e 


b =2r =3N. 
N o 


The approximation of V N N enables a solution of (3.7) to be 
analytically determined: 


( 0 ) 

< V r > = 



(j) J,(S7r) r<r, 

[A AaA -l 


where __ 

<V = 

for k =© corresponding to the usual procedure of multiplying the 
solution by a k^ dependent normalization factor to properly take 
care of both the k N dependence and the normalization to a plane 
wave at infinity (see below) . 


On matching function and derivative at r=r o and using the 

8 

well known asymptotic form of , one obtains to leading order 


, (<f ) \ ^ 
U U (r) = 


r 


>/4 


C<KL' ( ^ rF) 


r< r c 




7a 7 


vjn (7- r u ) A 


cv* (<*^ (r , ) f 3 -f) 


(3.9) 

(3.10) 

( 3 // 4 ) 
r> r 0 

C 3 , !2) 


(for r> >t ) 

( 3 . U b) 

r>r 0 
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IV. APPROXIMATIONS OF THE OPTICAL POTENTIAL 

We consider here three approximations of the optical potential, 


(i) The first includes only the lowest energy term coming from 


the 


<P — Q $ — <p (r t ) <p (tl) 


~Al+i 1 'N-h A7-/ 

Substitution of this into (2.15) gives rise to an integro-dif feren- 

tial equation. 

f d X _ V (r) 1 ulJ r) - R (r) V (r) < w |/ . 

[ rfyx, A//[/ J ” V+l W+) 


{4.1) 




't-/ r * — 0 


In this case because we have a separable kernel , the solution is 
given by 

U A/ fr )~ U J <J ' (r ) f C U J” fr ) ( 4 ' 3 ) 


where (r) is the homogeneous solution Eq. (3.11), is a 

solution of the homogeneous eq. 


fd Z __ 1/ (r) 1 LC Jr) - - // rr) /p Jr) 
l frx- Wa/ J J A hjtf+l * H 


and C can be solved for to be 

(o) 


C= -K 


E - «« + K 


TT)' 


(.1 #) 


(4.5) 


with 


and 


// fr) u ,(r)dr 



K ,n = f R fr) 1 / fr] u '%) dr 

J A'f-, A / /f f+, 


d,D 


H-7) 
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The coupling potential V N N + 1 , Eq. 3.6) is also a complicated 
function which can simply be approximated: 

y ft) £ o.g [/ +- Xfr/a) ']"*- (4. S’) 

2 

In Figure 1 we plot N V N N+ ^ vs. r for two values of N exactly cal- 
culated from (3.6) together with the approximation (4.8). 

- ;> The convergence as a function of N can be appreciated by our 

pointing out that the difference between N=10 and N=ll results would 
be indistinguishable on the graph. The fit of (4.8) is not perfect 
around r/r o =~.5, however, our results below are not affected. A 
better fit can be obtained with 

V = /3/V"7 Z /TJJTTr^T +( r /rcf + /*/£ 

hW+t L 



The solution of the equation, (4^4), is effected with a 

Green's function technique 




,'f) K fr>) Ur> 

L V// h./V-H 


where the Green's function is 


, (o) (o) 

G(r,r ) = (-2II)U^ (rjqr ( r> ) , 


( 4 . ■?) 


(4.101 


and is (any) irregular solution of the homogeneous equation. 

We choose the complementary form of (3.11) whose asymptotic form is 
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The details of the quadrature involved in (4.8) are given in 
Appendix A. The result is (r<r Q ) 

0^ 1) =U^ o) (r)II(r)+^ (o) (r)I(r) 
where I (r) and II (r) are given in (A. 7) and (A. 10) . 


The the N dependence of is exceedingly 

simple. One finds 

K<0) - iU 

The N dependence of is derived in Appendix B: 


K 


0 ) 


C j -f C ■) /lAsy\ { 2 - i/ZlP ^ . 7T~) 


A 


,/■ 2 — 


In order finally to evaluate C of Eq. (4.5) and thus 

U^j(r) of (4.2) we need to know the energy differences E- . 

The total energy, as was indicated, is taken as that energy to 
th. 

excite the N level 


E-E = " ~ 2 ' 
KT 


and from Table I we find that can be well fit by 


£ 


N+l 


1.27 
(N+l) ‘ 


To lowest order therefore 


>w 4 + i ■ 41 + 0 <1 > 


The function (r) in the region r~r Q is dominated by the term 


(r)I(r) by noting that for r~r Q 


I (r) « £. 

N 


i/2 


(4.12) 


(4.13) 

&■!+) 


.15) 


4.16) 


4.17 


.18) 
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as opposed to II (r) « N [using (A. 11) and (A. 12)]. Thus 

putting these behaviors together we find 


jtL- 


A/ — > oo 


u A / r y 

/V 


f, <«) 




where 

^(N) = CN 

l-Bsin(2l^l2 7 Nir) 




(4.20) 


The above is the essence of the P\l/ contribution to the wave 

■N 

function, however the total wave function includes a contribution 

QiJj . This may be derived from using the relation 
-N N 


= 1 QHPtf) 

e-qhq 


(4.21) 


(5) 

Eq. (4.21) is the first step in deriving the optical potential' 
equation (2.9) from the Schrodinger equation (2.1). For the one 
term approximation that we are here considering, (4.21) reduces 
to 





(4.2 ft) 


where is given in 


. 1). The integral reduces to 


<?> /— / [ r , V. _ /Of fr\ \/ (r) OL^fr)^ (4-2 2) 

/Vf/ < / ~p= J e *r+t w 


and using (4.19) for tl^r) reduces this to a form involving K 

(1) 


(o) 


and K 


One finds in fact 


k/bere. 


9 V) - j, /A/> d^ +/ %) 

%, Wj = AJ Y K s/TT p fa) Ar °'J ft. Z3 ) 
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which upon substitution reduces to 


a (fif) - 4/ ^ f C] ~ Q. Th j (4-2-4) 

where the C 5 and B '5 are constants which can in principle be de- 
termined. 

The threshold law is derived for this as wfill other approxi- 
mations of the optical potential in the next section. 

(ii) The second approximation we shall consider is motiva- 
ted by the observation the optical potential (2.15) is (formally) 
dominated by states = E (dominant pole approximation). The 
actual states for which 

= E n (4.25) 

are readily deduced from (4.15) and (4.16) to be 

11 ■ .12 N (4.26) 

[Cf. below Eq. (2.17)]. In other words we consider an optical 
potential based on one intermediate state 

< 3 ? = -L j? (r,) fa , 

2 

Because the energy denominator vanishes (to order N ) in this 
case, p(N) may be simply gotten by putting E= <f in (4.5). Then 
using (4.19) we see that 

/ (9) 

< — &■***) 

C, / Q .^ /2^/t/7r; 
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On the other hand q 2 (N) must be evaluated more carefully, 
because the expression (4.21) is indeterminate. One finds 




C 


K_ 


fo} 




N 


3 A 


C. f -h (iJJi A/) 

The relation (4.19) for U^(r) applies in this case also. 

(iii) Finally we consider a closure approximation; the in- 
termediate energies in (2.15) are replaced by a mean energy 
so that 

n j ' Q 2 

^ - 1_ 2 ? HQ* Nv ><Q* N v HP = 1__ PHQ HP 

E -C e-c, 

N V 


= 1 




_ PH (Q) HP 


TV 


* 1 PV(1-P)VP ^ 

E “^v 

since [P^H^I = 0 - PQ. 

In the uncoupled approximation, P reduces to 
p -> 

The optical potential of (2.10) becomes 


1 

op 


E -£ 


[PV 2 P-(PVP) 2 ] , 


TV 


where 


V= - 2- + 2 


r / r, v- r. 




(4.29) 


(4.30) 


(4.31) 
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Based on an approximation similar to that used to derive the 


form (3.7) for V ^ , we can show (Appendix C) 




V r / 


n*?A> * 




^>r I; 


(l 33 -) 


Therefore, with use of (3.7) for V NN , (4.30) becomes 


9 




N 


a 

4- 


£-f r.4 

_ s v i 


r,<r 0 


*^/>C 


( 4 < 3 3 ) 


The N-dependence of may be estimated variationally (Appendix 
C) to give 

£ a i 
N ¥ 


(4.34) 


Our approximations are not sufficiently accurate for either the 
sign or the N-dependence of the energy difference E N - ^ occurring 
in (4.33) to be determined. The best we can do is to limit the 
difference by 




(4.35) 


This gives rise to an effective local potential for the scattering 
function 


a? + Vc<r>1 Vr) = ° 

where 



(4.36) 


3 7) 
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In (4.37) we have made the assumption that the optical potential 
in fact exceeds the dipole potential in some finite region (i.e. 


r^>r Q ) , where r^ may be determined by the condition 




which leads to 




If (4.39) does not lead to r Q >r , then the equation and solution 

p o 

revert back to of Eg- (3.11b). 


The solutions of (4.36) must again be determined by matching 
and one finds to 'lowest order r<r 

o 

(r) - Ug, (o) Xr) (4.40a) 

and for r Q <r<i^\ 



P 2 = 6^ (4.41) 

iv5i 


The solutions for r>r^ go into the general form of the dipole 
potential given in the lower part of Eq. (3.11a). The coefficients 
are again determined by matching / but we shall not consider them 
further. 
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V. THRESHOLD LAWS 


Threshold laws are calculated from the expression 

^ t is the yield as a function of the available energy E after 
ionization. M.E. is a matrix element 

M.E. = < X IV I 

which causes the trans 2 tion from the unperturbed initial state 



(rz) 


$ t =- A<r ( /Ar, ) ( r i) (tt) 

to a final state \b. y the calculation of which we have discussed in 

-hSF 

the previous sections. The quantity gis a normalization constant 

<r- ~ 

which adjusts the 11^ to be a plane wave at r^-*- 00 . rg was evaluated 
in the appendix of Ref. 1 

-y = r 

/ / /'■***/- YVVWoV) -fry Cr. *) 

It should be noted the factor was omitted in Ref. 1. 

N 

9 

(We are indebted to A.K. Bhatia for finding the error .) From 
(3.12) we see that 

iLn c/ ~ jy a/ (sir) 

A/- > / Y 


and in (5.4) 




where r m a matching radius beyond which only the dipole potential 
2 2 

b^/r and the outgoing energy k N enter the equation for U^. The 
point is that the k N dependence of is completely absorbed in n 
and the calculation for are done at k N =0. 
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In order to arrive at the ionization region we have assumed 
an analytic continuation of the inner electron's energy from E^= 
[Eq. (2.7)] to /v =+l 


N 


N 


E N' >Ar N“^T 

N 


(5.7) 


This continuation is motivated by the well-known fact that a 
Coulomb wave for negative energy becomes a positive energy solution 

g 

by changing N-+— i/x*J^, in the confluent hypergeometric function . 


The threshold laws are then derived from (5.1) wherein from 
(5.4) the explicit k N dependence cancels out, and the remaining part 
of the integrand is converted to a function of via (5.7), so 
that integration gives the E dependence of fa, which is what we are 
seeking. 


To gain confidence in the analytic continuation - let us 
consider as an example, the homogeneous solution of the Hpp prob- 
lem, i.e. given in Eq. (3.11b). Here the matching radius is 

naturally taken as r Qo 

r = r ( 

m o 

so that to leading order 

( c) 

(X Uy , & ' f/fr Q - yrrfa) 

dr £ ^ /.? jrfj ) V* ' 

and 

)- /L tcuyx 4 ,r ¥n - ) 
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Thus using (5.5) and (5.8) for the N-dependence of r Q and a N we 
find 



c - ^rj c^- 

or finally, 

1 

What is nice is that the oscillating factors in the denominator 
cancel away; we shall find this to be essentially always the case 
as regards r| . 




'A 


>% ^ A 







The remaining piece of the integrand is the matrix element 

which in this approximation is 

!■: ) 


M °- < I vi 


In Appendix D we show 

M on N -3/2 
o 

With M.E. being M q in this case and substituting for n , we find 


or 



fs:?(b) 

(5.9b) 


X Q WE (5.10) 

A linear law is precisely what we expect in this approximation in 

which the potential felt by the outer electron is purely Coulombic 

on the inside and attractive dipole on the outside. For it is now 

well known that the latter also causes a finite inelastic cross 

to 

section at threshold y and this is guaranteed in our formulation 
by the normalization constant n . (The subscript on will attempt 
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to specify the particular approximation used.) 


We next consider the lowest term and dominant pole approxi- 
mations of the optical potential. In these cases the matrix ele- 
ment contains a part from Q-space (the term multiplied by q^(N) 
below) in addition to the P-space contribution: 

M.E . = El+N _3//2 q i (N)jM Q +p i (N) . (5 . 11) 


The index i=l,2 specifies the two approximations. is the part 

of matrix element coming from the irregular solution part of tl^: 

M 1 =<U fj (1) (r.H (r 2 ) |V|4 > (5.12) 

initial 

r l 

(r) is given in (4.12). Although is more difficult to 

calculate exactly we have shown in Appendix D that 

M < i (5.1 Zb) 

In 


a-mcA- (M) 

From (4.20) and (4.28a) we see that p.^ (N)^ are essentially 

proportional to N thus p. (N)M. is smaller than the M term of 

i l o 

(5.11) . Concerning the evaluation of ti the 

dominant term of is dominated by p^(N) [Eq. 4.19] which one 

power N larger than in the case.' On the other hand the 

logarithmic derivative is the same 



aside from the interchange of sine and cosine factors. The same 
l yiiercbA-* (f ) , . 

A is true v for vs. j therefore the oscillating factors con- 

tinue to cancel out and we are left with 
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*1 00 


(5.14) 



In comparison with n this normalization constant is dominated by 

the N in the denominator which causes the threshold to be contain 

/A 

an extra power of E' 

^“/Ml-Bsin (7-f^r) ] 2 dw 

which to leading order is 

(5.15) 

This result is at first sight very unexpected. However, from 
the point of view of the lowest optical potential term approxima- 
tion, wherein we have shown that this term is rigorously repulsive, 
the result is seen to be reasonable consequence of the repulsive 
optical potential term retained. In the dominant pole approxima- 
tion, in which the term selected is at the border-line between 

•w 

attraction and repulsion, the physical origin of the result is not 
clear. This is particularly true because the shift, of 

Eq. (4.14), is also very likely to be an oscillating function of 
N. [We have, together with Dr. Bhatia, numerical solutions of the 
exact lowest term eqs. up to N=9 which indicates that this is the 
case.] The lesson to be learned is whether we understand or not 
the optical potential can be expected to have a profound effect 
on the threshold law. 

We finally consider the closure approximations. Here we have 
the. possibility of many results in view of our ignorance of the 
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sign and the exact N-dependence of E^- even within the confines 

of (4.35). We shall subdivide these into attractive and repulsive 

2 

cases, both with the assumption that the 3 / r 4 potential is 

2 

stronger than the b. T /r in the region r <r<r„ [i.e. r„>r from 

IN Op p O 

(4.31)]. 


In the attractive case we find, using the upper solution of 


(4.40) , 



(S. /£>) 


Furthermore we have shown in Appendix D that the N dependence 
of the matrix element is not altered by the contribution of 
U^ 6) from r Q <r<r g providing E N -^eC-N -2 . 

M.E.°c M c , oC //"V 2 - (5.17) 


Thus substituting gives 




(5.18) 


(5.19) 


which implies from (4.39) 

r|p*J 2 , (5.20) 

gives using (5.7) in (5.18) : 

3-L^E (5.21) 

Another conceivable alternative would be for example, 

— ^ 

E n - • For this case the matrix element would be dominated by 
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the Q¥ part of ¥, as shown in Appendix D. The net effect would 
be to give an E ly/4 threshold which we shall not pursue further. 


Penultimately we consider the repulsive closure approxima- 

/ O \ 

tions corresponding to v ' of (4.40b). Here the normalization 
constant turns out to be 


')? A 7 

'% f j3 ) I f 

The cosh factor in the denominator which appears to dominate n 2 is 
however, cancelled by a similar factor in the transition matrix 
element (Appendix D) 



ne. * Co z<^_ A, r/s 

8 Y 

Using -^CN 1 , where y>o in all cases, we are left with 
o 

£ 


j. 


oU 1 


u 


■x ^0 vp 

If now we restrict ourselves to quadratic dependence of r^ on 
N specified by (5.20) [albeit now in a repulsive sense], we find 
X at E 5/2 (lnE) 2 


(5.25) 


There is absolutely nothing at this time which prevents the 

effective optical potential, as contained in the energy difference 

E N~ “tj from being an oscillating function of N in sign. From (4.39) We 

/(tKa^ r when E N ~ changes signs. Assuming the amplitude of 

2 

these oscillations is 1/N , we see from (5.18) that attractive 
portions give a linear rise whereas from (5.24) the repulsive 
portions are essentially flat. There are of course many analytical 
functions which can give this type of behavior - an example would be 
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(5.2 6) 


ce [1-Csin (InE) ] 

where in order for the slope not to be negative we must have 
C<2 -1 / 2 . 

A sketch of such a threshold law is given in Figure 2 for C=l/2. 
It can be seen that such a threshold is distinctly non-linear. 

In addition its oscillations about the line ££.=E continue right 
down to origin. 
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V. DISCUSSION 


We have not attempted to derive a unique threshold law. 

Our purpose in this paper has been to present what we believe is 
potentially useful and rather different approach to the problem. 
The approach naturally leads to the optical potential as the key 
element beyond the obvious potentials that the outermost 
(scattered) electron sees. We have been able to show rigorously 
that this optical potential contains repulsive terms, although we 
have not been able to determine whether the repulsion or attracti< 
dominates in the potential as a whole. The repulsive approxima- 
tions can lead to a considerable diminution even beyond a simple 
phase space dependence on E : 

^ vcfdtmK fc-ti'-io * / 

/Vut oul J ' > is 

Conventional wisdom on the subject might have dictated that we 

delete those approximations which lead to a higher power than two 

other 

however, we have included them because we know in, contexts that 

A 

threshold barriers can have^overwhelming effect on threshold 
cross sections and we cannot exclude that situation here. 

, 11 12 
We have not discussed the salient recent work ' which 

attempts to justify the Wannier law on the basis of a more con- 
sistent WKB approach. That work is significant but it is not 

rigorous. It can be shown in fact that the Wannier threshold 

4 -1 

law remains the same m that approach when the r ^ interaction 

with that of our model (r^t^) ^ • However, the most provocative 
of our results is the oscillating derivative threshold expressed 
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in (5.27). It is the possible existence of such a threshold law 
which makes a reliable calculation an attractive initial endeavor 
as part of the general problem of synthesing the optical potential 
in a definitive manner. 
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APPENDIX A: Evaluation of 


U) 


We wish to compute the function ^ of (4.9) with 
G(r,r') given by (4.10), Uj J (0) by (3.11), v N (0) by (4.11), 

V N N+ ^ by (4.8), and by a similar asymptotic expansion? 



C N is the Nth (last) coefficient in the expansion of the 

g 

Rjy (which is rxR^ in the notation of Bethe & Salpeter ) 



(-1) N - 1 2 N 


N 3/2 N! N 1 ^ 1 


(A.2) 


It should first be noted that our approximation of R^r) 
is not continuous at r=2r Q and that the part for r>2r Q is 
the very asumptotic form to be used only in showing that 
contributions to (r) from r>2r^ are negligible. (Cf. F r i). 


It is to be emphasized that the rhs of (A.l) is divided 
into two regions at r=2r Q . The fit of by the rhs of (A.l) 
is no longer accurate even for r=r Q ; the well-known but com- 

g 

plicated WKB expressions for could be used between the 

4 2 

classical turning points 0, jr Q =2N , particularly around 
r<rg and they are much more accurate that (A.l). However 
reference to Figure 3 shows that the exact function is some- 
what larger than (A.l) around r=rg and oscillates more 
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slowly, it does continue to oscillate beyond r=r Q but it 
ceases to oscillate and is much smaller than the rhs of 


(A. 1) at r=2r Q . For this reason we believe, for integra- 
tion purposes, these compensating effects are adequately 
accounted by simply continuing the rhs to r=2r Q . 

From (4.9) 

“v* =■ ' 0 >rr) Jrr> * u J 0 l H fr ) (A. 3 ) 

with 

^ r 

Kr) = inf uf'(x) Rjf) [/ m cJ . X (A.fJ 

and J ° 1 

TL( r) - xrf rf'n) tjx) KJf ^ (4 s) 

yc 7 

Assuming r<rg and using the equations stated above, we find 

r'f 

2(f) C <rv' (j?* - m ) d* CA-i) 

4 Cr^ixt 7p ’ 

Replacing cosine square factor by its average value (1/2) 
gives 

T(r) * />/ ■+- jirf e t(Ln '/% j 0-7) 

For r<r^ the factor II (r) contains two contributions. 

We shall show later that the contribution from 2r^ to 00 is 
negligible, therefore we have contributions from r to rg 
and r Q to 2r Q . 

The first is 
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Making similar approximations as above but retaining the 
sinusoidal factors, we get 


7 Z(r) ~ 

< 


(Si f 




The factor jis a number of the order l^l^lO to make up for 
the fact that the bound 

r^ 2 - > (<~o *-27-) f 1 ' 


fA.j) 


has also been used in deriving (A. 9). Here the sinusoidal 
factor cannot be dropped because its mean value is zero, but 
(A. 9) can be integrated to give 







The expressions (A. 10) and (A. 7) into (A. 3) are to be 
used in Eq. (4.12). Note that the N dependence of II(r) 
dominated by the second term and that (for r^rg) it is 

iKrlotV/, (AH) 


We shall now show that the contribution to II (r) from 

-3/2 

rQ<r<2rQ is of maximum order N ' . Using the r>rQ 

form of v~ N (0) gives for Iljr) ~ 

* r o 

£^/U^lx v ^ _3Tr J ^jdr- I 

0r<+zr)'~ / 






(C) 


'C, 
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.IT . 1 (ir) 

The potential 0 being bounded by — and R* T1 , , being 

2 r o N+ ' 

(r 0 +2r)^ U 

bounded by setting the cosine factor equal to one: 

< far) '^ / / jr/v^ 

A'-hi / 


gives aside from numerical factors 

) y 


/ R V r 

1/ ^ A/./Zbt 

• o / 


$1 1- , 




^ /? i Cru^jt 


where. we h^Ve l<?£ . h/e ca 77 

tn-teyr*/ cnthe rhg h <^> J since the lay < 
co nfrL l u htc-ry corves f ee m ^ j ,„*//. Th us 


wJ C) P 1/ Jr 

» a /A/ 


-c -~_ 


KS7 -i [ PS t /tA^ / <w_ 


^ ^ ? T ~ 

7 " h t S ts /V de.J>e. y\ c4e y\ c e 4 S 

u/ <? ^ k ti / e-f h tv' l irT 

7 Ur) - Hjr ) = Chih-) 
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We next show that the contribution to II (r) (and also for the 


similar contribution to I(r)) from 2r o <X <c0 is truly negligible, 


. oo 




(A. 1 ) , 

CO 

Vw 

( 3 .ii) ©wy.//j iveyv-M' 

0 . (c) 

U-H - 

V- 

\ d f * Q+, 



Vi 


1 






X - r /W7 dr 

r~*~ 


Use of r a N z shows the rhs is bounded by 
o 


o* 


£ 


/ 


e' A/*t ar 






Below we drop factors N v where v is any number independent of N 




which using Stirling's formula is 


oce N(log6-2) =;e -.21/V 


or finally 




A’Av a 7 ^ 


K 


'V 


w 




That is smaller than any inverse power of N. 
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APPENDIX B: Evaluation of K 


(1) 


From (4.7), (4.12) and the fits to V N N+ ^ and 
main contribution comes from 


, the 


J — 1 U « ^ ^ ^ H t - 

o (f 0 +2r) z [ J 


Consider the first term in curly brackets; using (A. 12) 


N h J T(r)cjr 



^ A ~/f c r ' /2 ' C(M - 

which is to leading order is ( ^ 'h 7 -'") 


which is to leading order is 
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Now use Dwight‘S, Eg. (416.17), 


Cct± 0-J> 


1 a kTl / L +- c e*(° _ coi-xp , , , 

-jf— 1 m - 


/ X -0L ^ 




and realize that the main contribution comes from first term. 


Thus find that the above contribution is 


aC /~Y'c ^ ^ 

r ‘ 1 ~T(uff^-o 


/ 



(f 


C 


(zJTTA/r) 


This is the order as the first term but of oscillating sign. 
The sum is 


C,_ iC, 

A/^ 


37 



APPENDIX C: Closure Approximation and Evaluation of the 

Average Energy <^ N . 


The closure approximation is introduced in (4.29) to 
simplify the Q-part of the optical potential . We con- 

sider here in more detail the evaluation of the average 
energy E which appears as a parameter in that approximation , 
(i) Since V assumes the form 


V = 




2 . 


2 K,_ 




f, + 

'Vr , 1 


we simply set 


- W- 




z > - 


and, from (3.7) 


r, 

A << t 


-CO 


- 4 A; 2 


. 7 /t 


tL 


'c_ 

to 




~ /oN 


M 


r, 


(C. i) 

> , ' 'c 

fc W 

«; > c 


4/r, 


U v ^ 1 -Ufa = IJ? 


<7^ 


, 7 F q 


/7^ 


therefore, * may be approximated as 


(c.y) 


n>n 


o r 


<? 


'i 



0 


fc <y> d**- 0( -f t ) 


A/ 


V 


r,+ 


% r Y^r) ^f r >~ r >) 


*7 


7>' 




fc.V) 
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where % is the unit step function. 


The approximations involved in (Ci2) is essentially the 
same as that employed in the evaluation of V NN and (C.3), so 
that is the form (C.4) is consistent with the P-part of 
the problem treated in Sec. / / ( . 


(ii) The evaluation of the average <£" is carried out by a 
variational procedure developed earlier. That is, we have 
replaced G® in (2.10b) by : 




9 


-> 


C,' 


e~cS 


9 

a 


ty(tz-f-j-) (y £ 

On the other hand, we can introduce a separable form for 


(C.r) 


with a set of variational functions Q<|> and Q(J>, as 


1 ^ Q* (c.6j 

It is to be emphasized here that the final state wave function 

we are trying to calculate corresponds to the elastic scatter- 
th 

ing from the N excited state at total energy E. Thus in 
the analysis below we will eventually put the initial and 
final state wave functions equal to each other. Consider 
first, however, the somewhat more general case where 
the transition- element is given by 

> <'- 7 ) 

with * 


X ■- 9 w y ? 
(xvnj 


(c. y) 
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sr 7 

We require that both (C.5) and (C.6) give the same 

u • j ( ‘ 

That is, (0 

<% X %->s. <%fq/E-E) X) 

which gives then the connection between c and < Cp cp 0 
Substituting for G® as given by (C.6) allows (C.9) to be 
solved for (fT in the form: 




<%<f?X<r?,7. r > ' 


<>/ 7y> 


rX<M 


Thus far, the trial functions Q<J> and Q4> are left arbitrary, 
except the normalization (linear) parameter which was elimin- 
ated by writing G® in the normalization-independent form (C.6) 
Now, we choose these trial functions such that (C.10) assumes 
a simple form, i.e. let 

Gf p 

Substitution of (C.ll) into (C.10) immediately reduces 
to a form 

J = . <% QHp XX- jL 

D 

where, using (C.8), we can write 


(C.ll) 


Cc.n) 


3 =<^ (,!!(> % } = <?> PV(f HQVPH'G 

' (c, t3) 


0 - 1 ( > =<7 , u 7 
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We can explicity estimate the N-dependence of B and D 
for using the result of Sec. Ill for the case P'J? = 

^ (r) 4 > n (^ 2 ) Firstly, consider the constant D, which becomes 

(using Q=l-P) as in (C.4) 

/ cXJ r . z I 2_ ■ / 

D =/ 4 l u v fr,)J -(fy,) 

) 

/ CX} 

^ Jr, tuSHifr,-* 

Using (3.11b) for U^ 1 "^ (the part for r>r Q ) and replacing the 

cos2 ( v ( r M ) - i/Wc -3 7T//J 

by 1/2 , we get for the integral 

J Jr V ^ / r -Vr ~ f r 

Thus 

2) ^ /^/ r 4^1/ 

The evaluation of B is longer and somewhat more involved. 
We have 



8 K Py F PVQ fj <? VP~y f > 


(C.14) 


(C. /S') 
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v/h<2v~e 


«?(«)/*, 

H - -l, +i^tir 

p = 5 ^VX(SAJ ; Q~/ - P 

Thus using h 2 <P N (2) =E N <p' N (2) gives 
QHQ = Q (h 1 +h 2 +V) Q 

=Qh 1 +Qh 2 -E N Q+QVQ . (C. 16) 

so that 

4 

B=£ B'i\ 

1=1 ' " y 

where the four terms come directly from the substitution of 
(C.16) into (c. 15) . Consider first 

B l =< U N ( ° )<!, N (2)VQh l V ' Cr N (0) ( 1 > < * , N (2)> 



- Q' (C . 17 ) 

Here we have used the definition 

Vo > x) V v > 


One can also readily find that 



42 



where in our approximation 


0 


r^a 


Af 


_z 

n 


r, { t >r o 


Each of these terms may &e evaluated in a straightforward 
manner except for the first term of (C.18) . In that case 
we use our approximation for V (but we neglect the cusp) 
before differentiating to find 

<vj’n) %(■>.) viyujh pv) 


Ayi 


~at 


[c //o) 



The first term of (C.21) cancels with the second term of 

_3 

(C.18) . All the remaining terms are of order N . Thus 

B oc-, (C.22) 

N° 

so that combining that with (C.14) , we get finally 
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APPENDIX D: Evaluation of Transition Matrix Elements 


We want to find the N dependence of M.E. in the 
ya t /\i'ou '5 approximations we have used. The P^ part of M.E. 
is 




(Dl) 


ff d<r < ^ U A/^) fyri) 

d 0 I > 1 


oC 

The first term of V gives zero by orthogonality and since 
R , (r 2 r 2 , the r 2 coordinate is confined to be close to 

origin; we can very accurately explnd 
2 2 2ri 


(D.2) 


r 3- +r 2 


r a 2 


Thus 

( M £) 



The lower limit on the integral can be extended to 0 (rather 
than r 2 ) because the integrand converges at the origin. If 
now we divide the integral into two regions. 
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we note that the first term is cut off by the oscillations in 


sinkr (which are independent of N) . And because (0<y<l) , 

, y 

r \ ^sinusoidal function of r; 


U N (r) i 


r 0 3/4j 


the second integral always converges and is proportional to 

N~^ / /the second term in (D.4) is negligible compared to the 
first term. This is true whether u N (r) | r> . r is either the 

attractive u^^ of (4.40b) or simply u^°^ (r>r Q ) of (3.11b). 
Thus the N dependence of (ME) is controlled by the first term 

r jp 

of (D.4) and this in turn is determined by <N|r 2 |l> which is. 


trivially. 


<H ' r 2l 1>: “ ^3/2 . 


(D. 5) 


M q is a special case of (ME) ^ so that we have finally 


M q « (ME) - |N _3/2 (D . 6 ) 

attractive r ^ 

We must also consider the contribution from the Q'F of 
the wave function. In the closure approximation (4.21a) re- 
duces to 


(p = L <p V p A (d 

B «-£, 

Assuming 

Q'ty - 6jn <p/z.) (b. 

' — A/ J / Vf. , 

and using PV - ) %<->-/ y ? = ~ ; 

one can reduce (D.7) to 


60 > 


V A/ A) U<<! 

A A/l I 4 / 
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To calculate the Q-part of matrix element 


fME) n ' I (,y ' ^’ /v) 

we bound the r<r Q contribution by 



c&u) 

Note that as long as I E N ~ £ N \e*Lff that both the P'F and Q¥ 

contributions to M.E. have the same N dependence. However 

if | E N ~ (£" N | <0 (— 2 ) then the QV contribution dominates. 

1 1 N ' 6 

8 2 

We next consider the repulsive ^ case which is now 
dominated by the contribution from r to r D . Using (4.40c) 

Op 

in the second term of- (D.4) 
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Integration by parts give 


Y 




c ll c^l l Mr f u )] - 3nr cv&k [jlf r'~ r u ') 


r 


r „ 

i dr *l r WO/i"'- C V 
Jr r 1 / 


( P- 14) 


% 

-2 -1 

In the region r>r , r <<r and since sinh is less than 
cosh throughout the interval, the second integral has higher 
inverse power of N dependence than the first term, so that 
we obtain in leading order 



Finally, we consider the part of the matrix element 
coming from the which occurs in the lowest term and 

dominant pole approximations : 

-J J SUL'Txir, % t (H) C-Ur chl (l)./c) 

i ■{. ' 2 - 

where U^ (1) is given by (4.12). The functions I (r) and II(r) 
can be shown to be of the order of or bounded by 



i %n) 


c < r a 
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Thus the two contributions to are 


yj lt »c <A/M |> J'V- A/ k (r)± f r) 2. ,W/r 


oC JL 




,J fr) J~(\r) jL Ajuvx 4r dr 


^/ 2 _ 00 ^ TLfr) st^ '^r cJr. 

Considering the latter first and using (D.18) » 


/ «*} 

6i/^/ /d^w ^ 


The integrand^ bounded at the origin, since both u^ ' and 
singer,) vanish at r=o, and it is bounded at 00 since f j ~ f Z*', 
Therefore, 


Mn_« A/ 


For we have 


a / y* 


' **' (f JIM dr 


If we use N ' [1 - e 


- 1/2 


] to interpolate on I(r) from (D.17) 


and put all sinusoidal factors equal to 1 , we can bound by 


M < 

'll . 


— ( - (L 

v-1- L 
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The term in square brackets forces the contribution from the 
lower limit of the integral to be 0 ; thus the major contribu 
tion comes from the upper limit, so that we are left with 
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: Comparison of£ N (inRy) With Various Thresholds 


H 


a 






MODEL ( V =| - r — ) 
r l r 2 


S-Wave 

( V =| ) 

r 12 


N 

^ n" e n 

^ N _ E N -1 ^ N _ E N -2 

^ N - E N -3 


< f_ E 

N a N-l 

£ -E 

N ^ N -2 

2 

-.0643 



-.0041 



3 

-.0296 



-.0127 



4 

-.0168 



-.0096 



5 

-.0103 

>0 


-.0074 



6 

-.00754 



-.0056 

>0 


7 

-.00555 



-.0043 


8 

-.004255 



-.0034 



9 

-.00336 

-.00008 


-.0028 



10 

-.00273 

-.00038 >0 


-:0023 



11 

-.00225 

. — . 00052 


-.0019 

-.00022 


12 

-.00189 

-.00057 


-.0016 

-.ooo34 


13 

-.00161 

-.00059 


-.0014 

-.00040 

>0 

14 

-.00139 

-.00058 


-.0012 

. -.00042 

15 

-.00121 

-.00056 


-.0011 

-.00043 


16 

-.00107 

-.00053 


-.00097 

-.00045 


17 

-.00094 

-.00050 


-.00085 

-.00041 


18 

-.00084 

-.00047 -.00003 

>0 

-.00078 

-.00041 


19 

-.00076 

-.00044 -.00007 


-.00070 

-.00038 

-.00001 

20 

-.00068 

-.00042 -.00010 


-.00063 

-.00037 

-.00005 

21 

-.00062 

-.00039 -.00012 


-.00057 
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Figure Captions 


2 / uk curves 

Figure 1: N V. T vs. r/r . The lower, are from numerical 

evaluations for N=2,5,ll as indicated. The analyt- 
ical approximation is the top curve. 

Figure 2: vs. E from Eq. (5.26), curve 2, (C=l/2) . 

Note that the curve is monotonically increasing 
but it oscillates (infinitely rapidly as E-*-o) 
about 5=E (curve 1). 

Figure 3: The radial function f^ N (r) for N=15 (denoted by *) 

vs. the approximation (A.l). The abscissa is 

*7 2 ^2 7 

P-C\ rj so that the average value of W Z-V 

occurs at p=N. Beyond r=2 r£ the approxima- 

tion in (A.l) is denoted by k and is barely 
distinguishable from the exact curve. 
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